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Self-consistent diagrammatic approximations to the Anderson or Kondo impurity model, us- 
ing an exact pseudoparticle representation of the impurity states, are reviewed. We first discuss 
the infrared exponents of the pseudoparticle propagators as indicators of Fermi liquid behavior 
through their dependence on the impurity occupation and on magnetic field. Then we discuss 
the Non-Crossing Approximation (NCA), identifying its strengths, but also its fundamental 
shortcomings. Physical arguments as well as a perturbative renormalization group analysis 
suggest that an infinite parquet-type resummation of two-particle vertex diagrams, the Con- 
serving T-Matrix Approximation (CTMA) will cure the deficiencies of NCA. We review results 
on the pseudoparticle spectral functions, the spin susceptibility and the impurity electron spec- 
tral function, supporting that the CTMA provides qualitatively correct results, both in the 
high-temperature regime and in the strong coupling Fermi liquid regime at low temperatures. 
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1. Introduction 

Over the last 40 years the Kondo problem has become 
an archetypical model of correlated electrons. The dis- 
covery by Kondo 1 of logarithmic anomalies in the per- 
turbation series for the electrical resistivity calculated 
within the spin exchange model of a magnetic impurity 
in a metal gave but the first indication of the existence 
of complex many-body effects in this and similar models. 
The key ingredient of these models is the availability of 
local degrees of freedom at the impurity, e.g. the sub- 
states of a local spin. Coupling of the quantum impurity 
to the conduction electrons of the metallic host induces a 
local interaction between the conduction electrons, and 
generates delicate many-body resonance states at low en- 
ergy. 

Soon after Kondo's discovery it became clear that per- 
turbation theory in the antiferromagnetic exchange cou- 
pling constant J is sufficient only at high excitation en- 
ergies/temperatures T so that JN £n(D/T) <C 1, where 
N = is the conduction electron density of states at 
the Fermi level and D the half bandwidth (high-energy 
cutoff). By summing up the leading logarithmic terms to 
all orders of perturbation theory, using the perturbative 
renormalization group (RG) approach 2 it is found that 
a single dynamically generated energy scale, the Kondo 
temperature Tk ~ D exp — jw^j ^ s generated. The per- 
turbative RG holds as long as \tl{T/Tk) ~ 1. 

In the opposite limit of T <C Tk, the impurity spin is 
found to be screened by the conduction electron spins, 
and in the case of impurity spin S = \ coupled to a 
single band of conduction electrons, a spin singlet reso- 
nance state is formed. This has been established first be- 
yond doubt within the numerical renormalization group 
method (NRG) pioneered by K.G. Wilson. 3 Later an an- 
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alytical solution of the energy spectrum of the Kondo 
model was obtained with the aid of the Bethe ansatz 
(BA), 4,5 fully confirming the NRG results. 

While these highly successful exact solution methods 
have led to a virtually complete understanding of the 
single-impurity Kondo problem, it became clear in the 
course of the mid 1980s and 1990s that the Kondoesque 
enhancement of electron scattering at low energies plays 
a central role even in more complex situations, like in 
strongly correlated lattice electron systems as well as in 
mesoscopic devices. In the former systems, this is because 
the strong on-site repulsion in conjunction with a very 
short correlation length effectively induces local moment 
physics, as has been shown formally by the mapping of 
correlated lattice models onto quantum impurity models 
with a self-consistency condition by means of the Dynam- 
ical Mean Field Theory (DMFT). 6 In the latter systems, 
discrete, localized quantum degrees of freedom, coupled 
to a continuum of conduction electrons, are often formed 
because of the spatial confinement of interacting electron 
states, e.g. in quantum dots 8 and carbon nanotubes, 9 
nanoscale constrictions and charge traps, 10 or magnetic 
molecules, 11 leading to Kondo behavior in the electronic 
transport. All these findings have demonstrated that the 
Kondo effect is an ubiquitous phenomenon in interacting 
electron systems and have made it a central theme of 
condensed matter physics. 

To tackle the complex quantum impurity problems 
arising in such systems, it is desirable to develop, be- 
sides the exact solution methods, flexible, approximate, 
but systematic techniques which do not rely on special 
symmetry conditions or a relatively simple model struc- 
ture, like NRG, or on intcgrability conditions, like the 
BA, and which are still capable of describing the high 
energy and low energy sectors of the model as well as the 
crossover region around Tk- In many respects, advanced 
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perturbation theory methods are attractive here because 
of their flexibility. They do not require a simple conduc- 
tion electron density of states (CDOS), and thus may 
be employed as impurity solvers for the self-consistent 
quantum impurity models of DMFT, 6 where the CDOS 
may acquire considerable structure, e.g. a Mott-Hubbard 
gap and a narrow resonance in the gap. They may be ap- 
plied to impurity models with more complex structure, 
e.g. Anderson models with several levels and realistic 
electron-electron interactions. They may even be used to 
describe quantum frustrated systems like multi-channel 
Kondo or Anderson models. An additional advantage of 
a method based on resummation of perturbation theory 
is that it may be easily generalized to nonequilibrium 
systems. 13 ' 36 The latter is of interest in the context of 
electron transport through nanostructures like quantum 
dots or point contacts showing the Kondo effect. 

The perturbation theory for Kondo or Anderson mod- 
els is complicated by the fact that the spin operators do 
not obey canonical commutation rules or that the elec- 
tron dynamics on the Anderson impurity are effectively 
restricted to single occupancy by the strong Coulomb re- 
pulsion, respectively. This difficulty can be overcome effi- 
ciently by representing the quantum impurity degrees of 
freedom by canonical auxiliary particle fields, 14,15 ren- 
dering Wick's theorem valid. For a description of the 
strongly correlated Fermi liquid (FL) ground state an 
infinite resummation of the logarithmic terms of pertur- 
bation theory is clearly necessary, where, however, the 
correct selection of terms is crucial. As will be seen, this 
can be achieved systematically by means of conserving 
approximations, derived by functional derivative from a 
generating functional. 

In the present article we review the method of con- 
serving approximations for Kondo-like models in terms 
of auxiliary particle fields, with special emphasis on de- 
scribing, by means of a single approximation scheme, the 
perturbative regime at high energies as well as the strong 
coupling regime below Tk- We review the Non-Crossing 
Approximation (NCA) as the simplest of a hierarchy of 
conserving approximations and discuss its strengths and 
its fundamental failures. A detailed analysis will reveal 
the origin of the latter. This will lead to the Conserving 
T-Matrix Approximation (CTMA) , which will be demon- 
strated to remedy the shortcomings of the NCA both in 
the high-energy and in the low-energy strong coupling 
regime. 

2. Model and renormalized perturbation theory 
for constrained dynamics 

In many situations the Kondo physics is more clearly 
described in terms of the Anderson impurity model, 7 

H = E^ c L c ^ + E^- ctB / 2 )4^ (!) 

ka " 

+ V^2(ct a d a + h.c.) + Ud]d^d\d l 

kcr 

where ct and d\. are the creation operators for elec- 

ku 

trons with spin | (spin degree of freedom a = ±1) in 



a conduction band state k and in the impurity level 64, 
respectively. To be general, we have included a magnetic 
field B acting on the impurity spin with Zeeman energy 
B = gjJLsB (/ib and g are the Bohr magneton and the 
Landc factor, respectively). The electrons may hop from 
the conduction band onto and off the impurity with am- 
plitude V. A sufficiently large Coulomb interaction U 
at the impurity essentially prevents double occupancy of 
the impurity. Provided that the impurity level is suffi- 
ciently far below the Fermi energy Ep, its occupation 
number will be close to one, meaning that a spin S = \ 
is located at the impurity. The Kondo spin exchange in- 
teraction model follows from the Anderson model in the 
limit of nearly single electron occupancy, after projecting 
out the high energy sector. 16 

Since V/D <C 1 usually, it appears natural to employ 
perturbation theory in V. This perturbation theory is 
complicated by two facts: (a) the impurity is an inter- 
acting electron system, for which the powerful quantum 
field theoretical methods like Wick's theorem, Feynman 
diagrams and renormalization of propagators and ver- 
tices are not immediately available; (b) the perturbation 
theory in V is characterized by logarithmically diverging 
terms, like the ones in the exchange coupling J. 

2.1 Pseudoparticle representation 

Problem (a) may be circumvented by working with 
pseudoparticle representations for the impurity states 15 
(or equivalently a resolvent operator formalism), 

4 = /t6W/-„, (2) 

where /J is the fermionic creation operator for the singly 
occupied impurity state with spin a and b\ a) are the 
bosonic creation operators for the empty and doubly oc- 
cupied impurity state, respectively. Since the pseudopar- 
ticle representation necessarily enlarges the Hilbert space 
into unphysical regions, care has to be taken to project 
onto the physical subspace, defined by all many-body 
states with pseudoparticle number Q = J2 a ft fa + b^b + 
at a = f. This can be done in an elegant way for any 
expectation value of physical operators acting on the im- 
purity states by working in the grand canonical ensemble 
with respect to the conserved charge Q and simply taking 
the negative chemical potential A of the pseudoparticles 
to infinity 14, 15 at the end of the calculation, e.g. for the 
impurity or d-clcctron Green's function in the imaginary 
time domain, 

Since (...) denotes the (time ordered) grand canonical 
average with respect to Q, this procedure allows the use 
of the full machinery of quantum field theory, including 
Wick's theorem. It is worth pointing out that the repre- 
sentation Eq. (2) with the constraint Q = 1 enforced by 
Eq. (3) is exact. Physically observable quantities are nec- 
essarily given by two pseudoparticle (or higher) correla- 
tion functions, which in principle requires the calculation 
of both, self-energy and vertex corrections. 

Inspection of the terms of perturbation theory shows 
that each contour integral along the branch cut of a pseu- 
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doparticle propagator carries a fugacity factor e~^ A and, 
thus vanishes in the limit A — > oo. As a consequence, any 
bubble diagram consisting only of pseudoparticle prop- 
agators vanishes in the physical subspace Q = 1, and 
any conduction electron propagator G ca (uj) appearing 
within an impurity diagram is not rcnormalized by the 
hybridization. Only in an expectaion value of impurity 
operators like Eq. (3), one factor e _/3A is cancelled by a 
corresponding factor in (Qe~^ s+X ^^ 1 ^), leaving a sin- 
gle pseudoparticle bubble. Details of the evaluation of 
pseudoparticle diagrams can be found in Ref. 17 

2.2 Exact properties of the pseudoparticle propagators 
The pseudoparticle propagators G x (lu), x = fa, b, a, 
are given in terms of the respective sclfcnergies (u>) as 

G x (w) = [GSM^-E.Mr 1 , (4) 

where the unrenormalized propagators read, after pro- 
jection, G° frT (uj) = [uj - (a + l)B/2}-\ G a b {io) = [to + 
e d - B/2}- 1 , and G a {u) = [uj - B/2 - We now de- 

rive their exact low-energy behavior. Rewriting Eq. (3) 
in terms of pseudoparticle operators and using Wick's 
theorem, it is seen that in G x (u>) the average (. . . ) 
is confined to the Q = sector because of the limit 
A — > oo. This implies that the pseudoparticle propa- 
gators have only forward in time propagation. There- 
fore, these propagators are similar to the propagators of 
the X-ray threshold problem, i.e. their spectral functions 
A x (oj) = — -lmG x (u) + i0) show infrared threshold pow- 
erlaw behavior, 

A x {uj) ~ <d{uj)uj- a * , x = fa, b, a uj<€.T k . (5) 

The exact infrared exponents a x can be determined for 
the spin-screened FL case by the observation that in the 
X-ray problem they are related to the scattering phase 

shift 5^ via a x = 1 — (s^/tt^ and by employing 

the Friedel sum rule. 18, 19 It links the conduction electron 
number in channel a, attracted from infinity in order to 
screen the impurity, to the phase shift, An^ = sjfK 

(x) 

Anca in turn is defined as the difference between the 
average impurity occupation number n da in the station- 
ary limit (time t — > oo) and the impurity charge created 
by the respective pseudoparticle operator, /J, b^ , or a) 
at t = 0. Hence, we have for G/ CT : AnW = n da — 1, 
An[ f J = n da >, a' ^ a; for G b : An { c b J = n da , and for G a : 
= n da — 1, a = ±1. This implies the exact values 
for the infrared exponents, 20 valid in the low-frequency 
regime uj < Tk, 

a f<7 = n d -n 2 d /2 + An da {B)-An da {B) 2 /2{&) 

a b = l-n 2 d /2- An da (B) 2 /2 (7) 

a a = -l + 2n d -n 2 d /2- An da {B) 2 /2 , (8) 

where n d — J2a n d<? IS the total impurity occupation 
number and An da (B) = n da — n d - a the difference be- 
tween the a = 1 and a = — 1 occupation numbers in a 
magnetic field. The dependence of the exponents on n d 
and An d {B) is characteristic for FL behavior because of 
the use of the Friedel sum rule (i.e. potential scattering 
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Fig. 1. Diagrammatic representation of the Luttingcr- Ward func- 
tional generating the NCA for U — > oo and the corresponding 
pseudoparticle selfenergies S^ CT , T, b and the impurity electron 
Green's function G da . Solid, dashed, and wiggly lines represent 
the conduction electron, the pscudofcrmion and the slave bo- 
son propagators, respectively. Throughout this paper, all lines 
are understood as the renormalizcd propagators, unless stated 
otherwise. 

only) in their derivation. Therefore, whether or not the 
exponents are reproduced by a given approximation can 
serve as an indicator of its correctness, telling whether 
the approximation captures the spin-screened FL fixed 
point of the single-channel Anderson model. In the fol- 
lowing we will focus mostly on the case of an infinite 
Coulomb repulsion U, where the more than singly occu- 
pied states (a, a') do not contribute. 

2.3 Conserving approximations 

The solution to problem (b) requires the selection and 
summation of the essential terms of perturbation the- 
ory, in pseudoparticle representation. Precondition for 
the projection onto the physical subspace is the con- 
servation of the local charge Q, implied by the symme- 
try with respect to simultaneous U(l) gauge transfor- 
mations of the pseudoparticles. Any approximation con- 
serving the projection has to preserve gauge symmetry. 
This is achieved by constructing conserving approxima- 
tions from a generating Luttinger-Ward functional <I>. 21 
The pseudoparticle selfenergies Ti x are obtained as func- 
tional derivatives of $ with respect to the corresponding 
Green's functions G x , and are thus functionals of the 
dressed G x . In this way a closed set of nonlinear coupled 
integral equations for the G' x s is obtained for any choice 
of <f>. The choice of $ furthermore dictates the calcula- 
tion of the physical impurity electron Green's function 
G da (and other correlation functions), since, by defini- 
tion, iVpGdcr = lini^oo 5$/5G ct7 is the single-particle 
conduction electon T-matrix. 

3. Non-Crossing Approximation 

In view of the small parameter V/ D it appears reason- 
able to start with the lowest (2nd order in V) approxima- 
tion for $. Since in this approximation, which sums up 
infinitely many self-energy insertions, the corresponding 
Feynman diagrams do not have any crossings, it has been 
termed the "Non-Crossing Approximation" (NCA) (see 
Fig. 1). The NCA has been pioneered by Keiter and Kim- 
ball using the resolvent operator formalism 22 ' 23 and by 
Kuramoto, who first recognized the conserving nature of 
the NCA. 24 It took more than 10 years, before the NCA 
equations were numerically evaluated 25,26 and further- 
more solved analytically in the limit of low energies at 
T = 0. 27 Most of these works were concerned with the 
limit of infinitely strong Coulomb interaction U —> oo, 
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in which double occupancy of the impurity is strictly ex- 
cluded. In this case, the selfconsistent NCA equations in 
conjunction with Eq. (4) read (compare Fig. 1), 

E& c » = Tf ^ f(s)Al(-s)G b (u, + e) (9) 

sr» = rW -./( £ )A» CT ( £ )G /CT (c + £ )(io) 

G^(u,) = J dsc-^[G f ^ + s)A b (e) 

-A f(T (s)G b {e - u)] , (11) 

where to is understood as uj ± iO for the re- 
tarded/advanced functions, T = irNo\V\ 2 , and A° c(7 = 
i ImG^/iVo is the local conduction electron density of 
states per spin, normalized to the density of states at 
the Fermi level No, and /(e) = l/(cxp(/3e) + 1) is the 
Fermi distribution function. For U —> oo the NCA cap- 
tures correctly the Kondo energy scale, and it provides a 
qualitative description of the formation of the Kondo res- 
onance. However, it fails in a magnetic field even in the 
high temperature regime T ^> Tk, producing a spurious 
resonance in the impurity spectrum at u> = 0, in addi- 
tion to the two Zeeman-split Kondo peaks. The origin of 
this failure will be traced in the appendix by means of 
an RG analysis. The NCA fails furthermore at tempera- 
tures T <C Tk, where spurious infrared singular behav- 
ior in physical quantities appears, in contradiction to the 
expected FL behavior. 17 ' 26 ' 27 The infrared exponents of 
the auxiliary particle propagators come out independent 
of n d , af CA = l/(N+ 1), a b NCA = N/(N + 1), with N 
the spin degeneracy, again in contrast to the FL behavior 
Eqs. (6-8). The low-T failure of NCA is less pronounced 
in SXJ(N) symmetric models with N ^> 1. NCA becomes 
formally exact in the limit N — » oo, with deviations ap- 
pearing in 0(^2). 26,28 Note that the deviation of the 
approximate NCA values for the pseudoparticle expo- 
nents ctf,b is of order l/N, not 1/A 2 as one may have 
expected. 

At finite Coulomb interaction U, the exchange interac- 
tion J acquires contributions from both, virtual excita- 
tions to the empty and to the doubly occupied impurity 
states. A simple generalization of NCA to this case, i.e. 
adding the second order perturbation theory for the two 
processes, fails to capture the simultaneous contribution 
of both channels in each order of bare perturbation the- 
ory, and therefore leads to a possibly by orders of mag- 
nitude incorrect value of T K (in the worst case it is off 
by a factor of Tk/D). An infinite summation of vertex 
corrections (the symmetrized finite U NCA, or SUNCA) 
is necessary to cure this problem. 29 

The failure of the NCA at temperatures T <C Tk 
can be traced back to its insufficient inclusion of coher- 
ent multiple spin-flip processes. The latter are respon- 
sible for the formation of the Kondo resonance state. A 
qualitative improvement therefore requires to include the 
proper vertex corrections which account for the domi- 
nant spin and charge fluctuations. There is strong evi- 
dence that this is achieved by the so-called Conserving 
T-Matrix Approximation (CTMA). 19,30 Before we dis- 



cuss the CTMA in section 6, we will in the next two sec- 
tions present relatively easily tractable generalizations 
of the NCA to adapt problems with several local orbitals 
and with a large, but finite Coulomb repulsion U. 

4. Multi-orbital Anderson impurity systems 

In contrast to analytical and numerical exact solution 
methods, the conserving technique is generalizablc in a 
straight-forward way for impurity models with several 
local orbitals. Such problems arise, for instance, in clus- 
ter extensions of the DMFT, 33,34 but also in transition 
and rare earth metal impurity systems. To give an under- 
standing of the wealth of the low-energy spectra of such 
systems, we here treat, on the level of NCA, a Ce impu- 
rity embedded in a metallic host, for which experimen- 
tal photoemission spectra are available. 35 Although the 
compound CeCu2Si2 is a heavy fermion lattice system, 
for which one expects in the lattice-coherent low-T state 
a dispersion of the Kondo quasiparticle resonance with 
width of order of T K , the photoemission measures the 
momentum integrated spectral density (implying only a 
slight broadening of the observed resonance) , so that di- 
rect comparison with the results of the impurity model 
are possible. 

Ce has seven 4f orbitals, which are spin-orbit (SO) 
and crystal-field (CF) split and have an overall valence of 
close to 1 in CeCu2Si2- Hence, one can assume a large on- 
site Coulomb repulsion U — > 00 between all 4f orbitals, 
and the model Hamiltonian reads in auxiliary particle 
representation, 

H = E^ C l C k, + E^mfUm, (12) 
k<7 mCT 

+ J2 V ™( c l bt f™ + h.c.) , 

mker 

where tdm are the energies of the 4f orbitals and V m 
the hybridization matrix element of the local orbital m 
with the conduction band, m = 1, . . . , 7. The resulting 
generalized NCA equations follow as, 

£ fmm'a{oj) = T mm , J ^ f(e)A c „(-e)G b (u; + e) (13) 
E 6 (w) = r ™'™/ V /( £ Mca(£)G/w^ + e) (14) 

mm' a 

G d mm>a(u) = J de C~ 0e [Gf mm > a (u) + e) A b (s) 

-A fmmlr7 (e)G b (s -u)] , (15) 

where T mm , = TrN V£V m >, and G/ TOTO / CT , G dmm >o- are 
the matrix generalizations of Eq. (4) associated with the 
selfenergy matrix S/ mm / CT in local orbital space. The fact 
that there are several single-particle levels at high ener- 
gies, ef m , grouped together in SO and CF multiplcts, 
implies a rich structure in the low-energy spectrum as 
well, because fluctuations of an electron from one local 
orbital via the conduction band to another local orbital 
are possible in second order in V m . Since there may be a 
spin flip involved, these fluctuations induce Kondo-likc, 
logarithmic divergencies in perturbation theory, how- 
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energy below E r [meV] 

Fig. 2. Upper panel: Theoretical T dependence of the 4f spectral 
function of CeCu 2 Si 2 for T = 15 K, 47 K, 105 K, and 200 K. 
The inset shows the calculated spectrum at T = 15 K . Model 
parameters: e/i = —2.4 eV, CF splittings of the J = 5/2 sextet 
A C f =30 meV and 36 mcV, SO splitting A S o = 270 meV, 
hybridization V = 200 meV. Lower panel: Photocmission spec- 
tra for the same temperatures. The experimental photocmission 
spectra, divided by the Fermi distribution to gain access to the 
states above the Fermi energy. The inset shows the raw data on 
the same energy scale prior to the division. All spectra are nor- 
malized to the same intensity at 100 mcV and are offset for 
clarity. 

ever located at energies corresponding to the differences 
Ae m i = e/ m — e/i > between an excited and the ground 
state local level, as can be seen, e.g., by inserting the bare 
pseudofermion Green's function Gf mma in Eq. (14). As 
a result, there appear (in addition to the central Kondo 
resonance at the Fermi level), multiple many-body reso- 
nances at elevated energies Ae m \ given by the SO and CF 
splittings, as well as shadow peaks at negative energies 
— Ae m i, the so-called SO and CF satellites. The height 
of these resonances depends roughly logarithmically on 
T, characteristic of Kondo behavior. In addition, the SO 
and CF satellites are broadened by the lifetime associ- 
ated with the inelastic decay of the excited local orbitals. 
The shadow peaks often appear merely as shoulders in 
the spectrum, as they correspond to transitions from an 
only virtually occupied excited 4f orbital into the ground 
state orbital. This physics has been analyzed analytically 
in Ref. 36 and is qualitatively well captured by the NCA, 
as a perturbative expansion of Eqs. (13-15) shows. The 
results of the numerical evaluation of the NCA equa- 
tions for the impurity spectrum -tr m \hnGdo{u — i0)] 
are shown in Fig. 2 in comparison to the photoemission 
spectra of Ref. 35 The central Kondo Peak D and the CF 
satellites (peaks E and shoulders C) with CF splittings 



of 30 meV and 36 meV are clearly visible in the main 
panels. The inset of the upper panel displays the NCA 
spectrum on a larger energy scale, showing the SO satel- 
lites (B and F) with a SO splitting 360 meV as well as 
the single-particle resonances A, whose SO and CF split- 
tings arc not resolved due to the large lifetime broaden- 
ings T mm i. The T dependence agrees qualitatively well 
with the photoemission spectra in the experimentally ac- 
cessible temperature range of T ~T K . 

5. Finite-C/ Anderson impurities: SUNCA 

For many quantum impurity problems, including the 
DMFT treatment of correlated lattice electron systems 
with a Mott-Hubbard metal-insulator transition, the 
large, but finite on-site Coulomb repulsion U is essential. 
However, a naive generalization of the NCA for finite U, 
the generating functional comprised of the two diagrams 
in the first line of Fig. 3, does not even give the correct 
Kondo scale for this problem. The reason is that in this 
NCA the fluctuations into the empty ("light" bosons b) 
and into the doubly occupied ("heavy" bosons a) impu- 
rity state are not treated symmetrically. The latter is, 
however, essential to obtain the correct spin exchange 
coupling 



via a Schrieffer- Wolff transformation. 16 Thus, the asym- 
metric treatment of either of the two contributions to J 
leads to a Tk which may be off by an exponentially large 
factor ~ e 7T ^ d ^ 2F or ~ e 7r l e d+ c/ l/ 2r j j. e . possibly by sev- 
eral orders of magnitude. For a correct treatment of both 
terms, there must be included, for each diagram with a 
light boson line, the corresponding diagram with the light 
boson replaced by a heavy boson line (which amounts to 
the exchange diagram of the former), and vice versa, on 
the level of bare perturbation theory. 29 ' 38 The importance 
of these vertex corrections has first been recognized by 
Sakai et al. 37 and later by Pruschke and Grewe, 39 with- 
out, however, formulating a conserving approximation. 
In these works a numerical evaluation has only been given 
in lowest order (the terms denoted UNCA in Fig. 3). On 




Fig. 3. Generating functional of the SUNCA. Zig-zag lines denote 
the propagators for doubly occupied or "heavy" bosons. The 
diagrams a) define the NCA for finite U, the sum of diagrams 
a) and b) the UNCA, which gives a symmetrical treatment of b 
and a lines to leading order only. 
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Fig. 4. d-electron spectral function for U = — 2e^. Solid lines: 
SUNCA results, 29 dashed lines: NRG results. 40 
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Fig. 5. Diagrammatic representation of the Luttinger-Ward func- 
tional generating CTMA. The terms with the conduction elec- 
tron lines running clockwise (labelled "spin fluctuations") gen- 
erate T^ c f\ while the terms with the conduction electron lines 
running counter-clockwise (labelled "charge fluctuations") gener- 
ate T( ci> ' (Fig. A-l, see text). The two-loop diagram is excluded, 
because it is not a skeleton. 



the level of renormalized perturbation theory, it means 
that for each dressed 6-line there must be included a lad- 
der vertex function with a-lines as rungs, and vice versa. 
The generating functional of the corresponding conserv- 
ing approximation, the so-called Symmetrized Finite-U 
NCA (SUNCA), is shown in Fig. 3. 29 The SUNCA is 
tractable with relatively moderate numerical effort, since 
it can be formulated in terms of no higher than 3-point 
vertex functions. The results of a fully selfconsistent eval- 
uation of the (i-electron spectral function within SUNCA 
are shown in Fig. 4 in comparison with NRG results. It 
is seen that the correct Kondo scale (width of the Kondo 
peak) is reproduced. However, the SUNCA solution still 
develops a spurious low-T singularity. This problem can 
be cured only by a more sophisticated approximation, 
the CTMA. 

6. Conserving T-Matrix Approximation 

To overcome the failures of the NCA to describe the 
FL strong coupling fixed point of the single-channel An- 
derson model, described in section 3, we have proposed 
early on the CTMA. 19 Before we present the explicit re- 
sults of the CTMA, we give the line of arguments that 
led to the construction of the generating functional of 
the CTMA. 



6.1 Construction of the CTMA 

As a minimal precondition to obtain a gauge symmet- 
ric description of the FL fixed point, a conserving ap- 
proximation must reproduce the correct pseudoparticle 
infrared exponents, Eqs. (6-8), whose dependence on n c i 
is characteristic for a FL ground state. It is easily seen by 
power counting arguments that any summation of a finite 
number of skeleton sclfcncrgy diagrams merely repro- 
duces the incorrect NCA exponents. 28 Hence, the gen- 
erating functional must be comprised of an infinite class 
of skeleton diagrams in order to describe FL behavior, 
in contrast to the Post-NCA considered by Anders (dia- 
grams up to 0(r 4 ) in Fig. 5). 31 Since the latter is a con- 
sequence of the singlet state formed at low T between the 
impurity and the conduction electron spins, one may ex- 
pect that higher than two-particle correlation functions 
need not be considered in the single-channel case. The 
approximations to the total vertex functions between 
conduction electons (c) and impurity degrees of freedom 
(pseudofermions /, slave bosons b) are then two-particle 
T-matrices. As the irreducible parts of these T-matrices 
we select the single (renormalized) b or / particle lines, 
since (1) in the Kondo regime these terms are the lead- 
ing contributions in the small parameter VNq] and (2), 
in the spirit of principal diagrams, this choice gives rise 
to the maximum number of spin and charge fluctuation 
processes in the T-matrices at any given order of (renor- 
malized) perturbation theory. This choice results in the 
ladder approximations T c ^, T cb for the total two-particle 
vertices shown in Fig. A-l (2), (3). The Luttinger-Ward 
functional generating these ladder vertex terms (and oth- 
ers) is shown in Fig. 5. It is comprised of all closed pseu- 
doparticle rings (skeletons) with each conduction elec- 
tron line spanning at most two hybridization vertices 
and has been termed the "Conserving T-Matrix Approx- 
imation" (CTMA). The pseudofermion and slave boson 
selfenergies, derived by functional derivative, are shown 
in Fig. 6. Note that the vertex equations for T c ^, T cb , 
coupled via the selfenergies, have parquet character. The 
analytical expressions for , E& and Gda are given ex- 
plicitly in Ref. 30 ' 32 

The selfconsistent CTMA equations are solved nu- 
merically by iteration. The results shown below support 
strongly that the CTMA correctly describes the strong 
coupling FL regime of the single-channel Anderson im- 





Fig. 6. The CTMA pseudofermion and slave boson selfenergies. 
In the second diagram of £ j a and the the third diagram of E(, 
the 1- and 2-rung contributions, and in the third diagram of 
£ j a and the the second diagram of £(, the 1-rung contribution 
to the respective T-matrices are omitted, (not shown explicitly), 
in order to avoid double counting. 
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Fig. 7. CTMA results (symbols with error bars) for the threshold 
exponents ctf and ab of Af a and At, for U — > oo, B = 0. Solid 
lines: exact values, Eqs. (6,7), dashed lines: NCA results. 



purity model. In the appendix evidence is provided that 
it also systematically describes the high-energy regime 
above Tk- 

6. 2 Results 

Pseudoparticle spectral functions. As a first indication 
for description of FL behavior within CTMA it has been 
checked if the CTMA reproduces the correct FL values of 
the pseudoparticle threshold exponents. 19 The exponents 
were obtained by a fit to the numerical low-T solutions 
for Af a , Ab in a log- log plot. The results are shown in 
Fig. 7, showing, within the error bars, good agreement 
with the exact values and especially the dependence on 
the impurity occupation number rid, characteristic for 
the FL fixed point. 

The static spin susceptibility of the impurity was calcu- 
lated from the spin dependent occupation numbers n da 
in a magnetic field B as, 



Xi(T) = 



dM 



(17) 



6=0 



where M 
and 



n da = lim -jr- — 

a^oo J dioe 



g^B^Ylia arLcr 1S t nc impurity magnetization 



iO) 



ImE, G Ja (u> - iO) + G b (uj - iO)} ' 

(18) 

XiiT) shows T-independcnt Pauli behavior for T ~ 
0.5Tr- and down to the lowest T considered, 32 indica- 
tive of the FL ground state with a completely screened 
local moment (Fig. 8). As expected, Xi(T) obeys scal- 
ing for at least a range of Tk within a factor 10, 41 
when plotted as a function of T/Tk and in units of 
(g» B ) 2 /(4T K ) = Xi (T = 0)/W, where W = T* K /T L is 
the universal Wilson number and T K is Wilson's orig- 
inal definition of Tk- 3 Since the BA and selfconsistent 
perturbation theory use somewhat different definitions 
of Tk, we rescale the latter in the CTMA as well as in 
the BA solution such that it coincides with T K . This al- 
lows for a quantitative comparison of B A and CTMA and 
for a determination of the CTMA approximation to W . 
The result is shown in Fig. 8, exhibiting remarkably good 
quantitative agreement. We obtain W { - CTMA ~ 0.462 as 
compared to the exact result W ~ 0.4128. 



a. 




Fig. 8. Static spin susceptibility as a function of temperature; 
BA, CTMA and NCA results (see text). Model parameters used: 
t d /D = -0.81, T/D = 0.2. 



FL behavior of the impurity electron spectral func- 
tion Ad a {u) and selfenergy Sd CT (w) is of prime inter- 
est especially for applications within DMFT. The imagi- 
nary part of the impurity electron interaction selfenergy 
S d(T (w - iO), calculated in CTMA 32 from 



(19) 



is shown in Fig. 9. T, da exhibits many features of FL be- 
havior. It has quadratic dependence on both, lu and T, 
at low ui, T, with no sign of a spurious low-energy singu- 
larity down to the lowest T considered (T ~ 0.01 Tk)- 
As discussed in detail in Ref., 32 the curvature of the 
quadratic behavior in u) and T is found to be in good 
agreement with the exact FL result, = a[ui 2 + 

(ttT) 2 ]/T k , where a is an exactly known prefactor. 32 ' 42 
However, the position uj of the minimum of ImS^u;) 




-0.2 



w/T K 



6 S 10 



Fig. 9. Imaginary part of the CTMA impurity electron selfenergy 
ImS SCT (w) for various T for the same parameters as in Fig. 8. The 
small solid arch represents the exact FL behavior at T = 0. 32 
The NCA result at T = 0.16T^ is shown for comparison. Inset: 
the T dependence of the minimum value of ImS SCT (black dots). 
The dashed-dotted line, drawn for comparison, has slope 2. 
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Fig. 10. d-clcctron sclfcnergy E<£ CT (u;) and spectral function 
Ad(u)) after the potential scattering correction has been applied 
(sec text). The arrow (lower panel) marks the unitarity sum rule. 



is incorrectly shifted to loq « —Tk, and ImT>da{u — i0) 
acquires negative values, thus violating the Friedel sum 
rule. When searching for the origin of this shortcom- 
ing, one must keep in mind that £d CT (u; w 0) is deter- 
mined via Eq. (19) by both ImGd CT (o;) and ReG^u;), 
and thus, through the Kramers-Kronig relation, by high- 
energy (potential scattering) contributions to Ad a (u)). 
Hence, the erroneous shift u>o may result from an un- 
precise calculation of Gda (w) at high energies, either nu- 
merically or due to neglect of high-order potential scat- 
tering terms. To correct this shortcoming, it has been 
suggested to add an appropriate, phenomenological real 
constant to Yida{u). Through sclfconsistency it acts like a 
chemical potential and shifts the minimum of IniE^o-^) 
to lo = 0. The results of this correction are shown in Fig. 
10. It is seen that by a single, real parameter, motivated 
by potential scattering contributions, the full FL behav- 
ior of E<j CT (a;) is recovered, and Ad(co) obeys the unitarity 
sum rule with good precision. 

7. Conclusions 

We have reviewed the method of conserving auxil- 
iary particle approximations for quantum impurities with 
strong onsite repulsion U. We have shown that on the 
level of generalized NCA relatively easily tractable ap- 
proximations are possible which give a qualitatively cor- 
rect description of, e.g., Anderson impurities with multi- 
ple orbitals or finite U (SUNCA) at not too low temper- 
ature, but invariably fail in the strong coupling regime 
below Tk as well as in a magnetic field. It was shown, 
both by numerical solutions and by a perturbative RG 
analysis, that the CTMA, although numerically demand- 
ing, provides a remedy for all of these failures, and gives 
an essentially correct description of the FL behavior in 
the Anderson impurity model. It remains to be inves- 
tigated if this makes the CTMA a suitable "impurity 
solver" within the DMFT approach to correlated lattice 
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Appendix: Perturbative RG analysis 

As seen in section 6, conserving approximations pro- 
vide a means for selecting the physically essential contri- 
butions even in the strong coupling regime, while pre- 
serving the symmetries of the problem. However, be- 
cause of the selfconsistent resummation it is non-trivial 
to find an approximation which also includes the lead- 
ing logarithmic terms at each order of bare perturba- 
tion theory. The latter is known to provide a quanti- 
tatively correct description of physical quantities in the 
weak and intermediate coupling regime, i.e. as long as 
ln(max[T, lo\/Tk) ~ 1. On the other hand, the leading 
log summation is equivalent to the perturbative coupling 
constant RG in the sense specified below. For finding the 
proper conserving approximation it is, therefore, useful 
to have a method which allows to analyze whether or not 
the approximation reproduces the correct weak coupling 
RG flow of the effective coupling constants under high- 
energy cut-off rescaling. In this appendix we present such 
a method and apply it to the NCA and to the CTMA for 
the Anderson model in the Kondo regime for U — > oo. 

The coupling constant RG starts from the observation 
that in a renormalizable system, like the Kondo problem, 
where the T, lo and B dependence of physical quantities \ 
is characterized by a single low-energy scale Tk (i.e. given 
by a universal function of the dimcnsionlcss variables 
T/Tk, u/Tk, B/Tk etc.), these quantities must be in- 
variant under high-energy cut-off rescaling D — > D — dD, 
and so must the two-particle vertex functions, the phys- 
ical quantities are comprised of. Hence, the direct calcu- 
lation of x is replaced by the renormalition of the vertex 
operators A of the Hamiltonian such as to keep the total 
(fully reducible) vertex functions 7 C ^ (spin and potential 
scattering) and j cb (potential scattering) invariant under 
cut-off rescaling, and subsequent calculation of \ t° lead- 
ing order bare perturbation theory in the renormalizcd 
couplings. The same RG scheme can be applied, if \ i s 
not given exactly, but in some conserving approximation. 
In this case, one must identify the approximate 7 C ^, 7 cfc 
as dictated by the approximation for \- Hence, the RG 
analysis of a conserving approximation amounts to the 



-4 







(2) 



Fig. 11. (1) Pscudofcrmion selfenergy expressed in terms of the 
exact 3-point vertex V (black triangle). (2) Representation of V 
in terms of the 4— point vertices 7 C ^ , j cb . The first diagram on the 
right-hand side is the bare 3-point vertex V. The pseudofcrmion- 
slave boson vertex 7^ 6 vanishes by the projection onto Q = 1. 
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following procedure, 45 see also: 46 

(1) Determination of the approximate vertex functions. 
In any approximation, a physical quantity can be ex- 
pressed in terms of the full hybridization vertex V and 
of bare Green's functions, the latter being RG invariant 
by definition, e.g. the / selfenergy (Fig. 11 (1)), 



s a-Ct s' 0)' -n 



/§/<«X* 



e + lo - iO) x 



(A-l) 



[V(«j, e + iO) G° ca (-e - »0) - ^(w, £ - iO) G° CT (-e + *0)] . 

V", in turn, is composed of four-point vertex functions as 
shown in Fig. 11 (2). This defines the approximate j c f , 
7 ch . Note that only for a conserving approximation this 
definition is unique, independent of the physical quantity 
considered. 

(2) Bare perturbation theory for j°f , r y cb . Due to the con- 
serving method, the four-point vertex functions are com- 
prised of the dressed propagators. In order to perform 
the perturbative RG, these must be written out in bare 
perturbation theory (non-skeleton diagrams). In leading 
log approximation, there are two significant simplifica- 
tions: (i) Each contribution to the pseudofermion self- 
energy contains at least one integration over the slave 
boson propagator Gb (compare, e.g. Eq. (9)). Its loga- 
rithmic contribution is cut off at high energies lo 3> Tk 
by the pole in Gb{Lo) at lo w ej. Note that this remains 
true even for the dressed Gb- Hence, only the bare pseud- 
ofermion Green's function appears in j°f , j cb . (ii) Selfen- 
ergy insertions to propagators which are integrated over 
in 7^, 7 ch are of subleading order and should, therefore, 
be neglected. 

(3) Perform on the terms obtained in this way the per- 
turbative cut-off rescaling scheme 2,42 to obtain the RG 
equations for the spin and potential scattering coupling 
constants. 

A.l Perturbative RG for the NCA 

Considering the NCA selfenergies in Fig. 1 (or the 
NCA d-electron Green's function) and following the steps 
(1), (2) above, 7^ is identified as the full slave boson 
propagator, ^ c J a , sty = \V\ 2 Gb, i.e. by the expression 
given diagrammatically by Fig. A-l (1). The conduc- 
tion electron-slave boson vertex function 7 cfc is by con- 
struction proportional to the unit matrix in spin space 
since the slave bosons don't carry spin, and can only 
give potential scattering contributions. Within NCA, it 
i s 7s ,b s = \V\ 2 G f,j5 ss '5 ss ' , i.e. of subleading order. Per- 
forming the cut-off rescaling D — > D — dD on 7^ , one 
arrives, to first order in dD and in leading logarithmic 
approximation, at the correction to the irreducible c-f 
vertex, 



ic/,(l) _ 



= -*o§ £[A (c/ Wv[A ( 



c/)l 



is'V.scr" 



with the spin indices as defined in Fig. A-2 (1). The dot 
in Fig. A-2 indicates that the integration over the con- 
duction electron energy is restricted to the infinitesimal 
range [-D, -D + dD], [D - dD, D\. 

A c ^, 8K c i are matrices in the product space of 
the pseudofermion and conduction electron spins, and, 




c c 



T(efl ! 
NCA , 




Fig. A-l. Diagrammatic representation of the Bethc-Salpctcr 
equations for (1) the NCA approximation to the fullly reducible 
conduction electron-pseudofermion p-h vertex "/^q A (2) the c-f 
T- matrix T°f , and (3) the c-b T-matrix T cb . The CTMA ap- 
proximation to f c f is the sum of (1) and (2). For the perturbative 
RG, all single-particle propagators are understood to be the bare 
ones (see text). The external lines are drawn for clarity and do 
not belong to the vertices. Single-particle reducible parts (SPRs) 
in (2) are omitted, since the SPRs of the total c-f vertex do not 
contribute to the RG flow because of rule (i) (see text). 



hence, can be decomposed into spin flip, Jj_ (S + <S> o~ + 
S~ <S> cr + ), spin z-component, J\\ (S z ® a z ), and poten- 
tial scattering, W (1 eg) 1), contributions, where S^, 
are the impurity and the conduction electron spin flip 
operators, and S z , a z the z-components of the impurity 
spin and the conduction electron Pauli matrices, respec- 
tively. Jj_, J\\, W are the respective running coupling 
constants, with J ± = J| = 2W = \V\ 2 G° h (Lu) = \V\ 2 /\e d \ 
for lo <C |e<j| for the Anderson model. This decomposition 
leads to the perturbative RG equations for the coupling 
constants within NCA, 

-JVo(J N J± + 2WJa_) (A-2) 



d\nD 



dJ\ 



II 



dlnD 
dW 
dlnD 



= -N Q (Jl + 2WJ ll ) 



w 2 



(A-3) 
(A-4) 



Eqs. (A-2-A-4) reveal two remarkable facts. (1) Within 
NCA the potential scattering amplitude W is erroneously 
renormalized under the RG flow, leading to a spurious 
divergence of W. (2) The RG equations are easily inte- 
grated to give (Jo is the bare coupling), 

1 + 2N Joln^ 

i.e. the spin coupling constant J as well as the potential 
scattering amplitude W diverge at the Kondo temper- 
ature T K = D e- 1 /( 2NoJo \ This demonstrates that for 
U — > 00 the NCA reproduces the Kondo scale correctly, 
however due to an accidential cooperation of spin and 
potential scattering terms. 

One may conjecture that this fact, that the NCA does 
not distinguish between potential and spin scattering, 
is the origin why the NCA gives a qualitatively incor- 
rect description of the Kondo resonance in a magnetic 
field, see section 3: While the divergent spin scattering 
amplitude depends on magnetic field and leads to the 
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describes the Zeeman-like splitting of the Kondo reso- 
nance in a magnetic field, since it is quantitatively cor- 
rect for the susceptibility as well (see section 6). 



Fig. A-2. Pcrturbative RG renormalizations of the (irreducible) 
c-f vertex. Diagram (1) is the NCA result, the sum of diagrams 
(1), (2), (3) is the result of CTMA, where (3) is the contribution 
from the c-b vertex, which is not scaling and is neglected. The 
black dots on the conduction electron lines indicate that the 
frequency integrals in those lines are restricted to the regions 
[-D,-D + dD], [D, D — dD]. 



correct Zeeman splitting of the Kondo resonance, 44 the 
incorrectly diverging potential scattering term is insensi- 
tive to the magnetic field and leads to the spurious third 
peak produced by the NCA at ui — at finite magnetic 
field. 

A. 2 Perturbative RG for the CTMA 

We now consider the coupling constant renormaliza- 
tion under the RG flow within CTMA. From the analysis 
of the selfenergies shown in Fig. 6 and of the correspond- 
ing Gda the CTMA approximation for 7 C ^ is given by 
the sum of the contributions shown in Fig. A-l (1), (2). 
The contribution Fig. A-l (2) leads to an additional c-f 
vertex renormalization under cutoff reduction, shown in 
Fig. A-2 (2), 

= +^of £ [A (C/) W,,UA (C/) U^ . 

s"cr" 

The contribution Fig. A-l (3) leads to the c-b vertex 
renormalization shown in Fig. A-2 (3), which, as a pure 
potential scattering term is not scaling, and, therefore, 
neglected. Adding dA c ^^ and dA c ^'( 2 ) and decomposing 
into J±, J|| and W contributions, one obtains the RG 
equations of CTMA, 

-2N J,\J ± (A-6) 



dlnD 



dlnD 
dW 
dlnD 



= -2N J 2 ± 







(A-7) 



(A- 



Inspection of the two terms in Fig. A-2 shows that they 
are the direct (2) and the exchange (1) scattering contri- 
butions and, upon collapsing the bare boson lines for lu <C 
\e d \ to the spin coupling, -\V\ 2 G° b {u) -> -\V\ 2 /e d = J, 
that they directly correspond to the two well-known con- 
tributions to the perturbative one-loop /3-function of the 
Kondo problem. 2,42 The failure of the NCA in the weak 
coupling regime is, thus, traced back to the fact that 
it incorporates only the exchange term, and the direct 
term is missing. This is cured by the CTMA: The RG 
equations (A-6-A-8) are identical to the well-known per- 
turbative RG equations of the original Kondo model. 
This proves that the CTMA incorporates the complete 
Kondo physics also in the weak and intermediate cou- 
pling regime, where the perturbative RG is valid. One 
may expect that, as a consequence, the CTMA correctly 
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